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ABSTRACT 

The hypernovae (HNe) associated with Gamma-ray Bursts (GRBs) may have a fairly 
steep energy-velocity distribution, i.e., E{> (3) ex /3~^ ior q < 2 and (3 > (3o, where /3 
is the velocity of the material and (3o ^ 0-^ is the velocity of the slowest ejecta of the 
HN explosion, both in units of the speed of light (c). The cosmic ray protons above the 
second knee but below the ankle may be accelerated by the HN shocks in the velocity 
range of /3 ~ /3o — 4/3o. When (3 < 4/3o, the radius of the shock front to the central 
engine is very large and the medium decelerating the HN outflow is very likely to be 
homogeneous. With this argument, we show that for q ~ 1.7, as inferred from the 
optical modelling of SN 20031w, the stochastic gyroresonant acceleration model can 
account for the spectrum change of high energy protons around the second knee. The 
self-magnetized shock acceleration model, however, yields a too much steep spectrum 
that is inconsistent with the observation unless, the medium surrounding the HN is a 
free wind holding up to a (unrealistic large) radius ~ 1 — 10 kpc or alternatively the 
particle acceleration mainly occurs in a narrow "dense" shell that terminates the free 
wind at a radius ~ 10^^ cm. 
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1 INTRODUCTION 

One of the most widely suggested sources of cosmic rays 
(CRs) is the supernova (SN) remnants (see Hillas 2005 
for a review). Since 2001, more and more researchers have 
noted that hypernovae (HNe), in particular those associ- 
ated with gamma-ray bursts (GRBs), may play an impor- 
tant role in interpreting the CR spectrum above its first 
knee\ i.e., ~ 3 x lO^'^ eV (Dermer 2001a, 2001b; Erlykin, 
Wibig and Wolfendale 2001; Sveshnikova 2003; Wick, Der- 
mer and Atoyan 2004; Wang et al. 2007; Budnik et al. 2008). 
This is reasonable since the average velocity and the total 
kinetic energy of HN outflows are much larger than that of 
normal SNe. So is the energy of accelerated particles. A re- 
hable interpretation of CR spectrum up to ~ 10^* eV thus 
should take into account the different energies and types 
among SNe (Sveshnikova 2003). 

How to accelerate protons up to an energy ~ 10^^ eV 
in the HN blast waves? Dermer (2001b) suggested that the 



gyroresonant stochastic acceleration might play such a role 
(see Fig. 10 of Dermer 2001a for a quantitative plot). Other 
authors (Wang et al. 2007; Budnik et al. 2008; Erlykin et al. 
2001) considered the self-magnetized acceleration model put 
forward by Bell & Lucek (2001), in which the magnetic field 
of the upstream has been significantly amplified by CRs. 
Considering the energy distribution of the HN outflow'^, 
Wang et al. (2007) and Budnik et al. (2008) suggested that 
with the second model the CR proton spectrum steepening 
around the second (first) knee could be reproduced. In this 
work, we point out one potential limit of such an interpreta- 
tion and show that the gyroresonant stochastic acceleration 
model does not suffer from that problem. 

This paper is arranged as follows. In section 2 we dis- 
cuss the energy- velocity distribution of HN outflows and the 
medium proflle surrounding the HN outflows. We find that 
for HN outfiows owning a fairly steep energy-velocity distri- 
bution, when Po < P < 4/3o that may play the main role in 
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^ The second knee and the ankle in the CR spectrum are at 



' 3 X 10^'^ eV and ~ 3 X 10^** eV, respectively. 



2 See Berezhko & Volk (2004) and Ptuskin & Zirakashvili (2005) 
for the influence of the energy distribution of normal SN outflows 
on the spectrum of accelerated CRs. 



accelerating the CR protons above the second knee but be- 
low the ankle (see the discussions below eq.(8) and eq.(16)), 
the radius of the shock front to the central engine is very 
large and the medium decelerating the HN outflow is likely 
to be homogeneous. In section 3 we calculate the change of 
the CR spectrum around the second knee which is caused 
by the energy- velocity distribution of HN outflows, and com- 
pare the results with the CR spectrum observation so as to 
constrain the models. In Section 4, we summarize our results 
with some discussions. 



2 THE HYPERNOVA OUTFLOW: 

ENERGY- VELOCITY DISTRIBUTION AND 
THE MEDIUM IT EXPANDS INTO 

2.1 Energy- velocity distribution of hypernova 
outflows 

HNe, especially those of which associated with GRBs/XRFs, 
are distinguished for the broad lines in their spectra, in- 
dicating very high expansion velocity of the ejecta. The 
modelling of optical light curves and spectra, in principle, 
can reconstruct the energy-velocity distribution of the out- 
flows. However, no reliable constraint can be given on the 
P > 0.3 outflow by optical even if that part had some opti- 
cal depth since the current optical modelling is not fully 
relativistic (J. S. Deng 2008, private communication). In 
SN 1998bw, SN 2003dh, SN 20031w and SN 2006aj, strong 
photospheric velocity evolution (the inner the outflow, the 
smaller the velocity) are evident (Hjorth et al. 2003; Maz- 
zali et al. 2006; Soderberg et al. 2006). The optical mod- 
elling of SN 20031w gave that material moving faster than 
(0.1, 0.2)c were ~ (1.4, O.1)M0, respectively (Mazzah et al. 
2006), implying a fairly steep initial kinetic energy distribu- 
tion E{> r/3) oc (r/3)-i -^ where T = (1 - (3^)-^^'^. But in 
other events, no result has been published. Soderberg et al. 
(2006) constrained the kinetic energy profile of HN outflows 
in a more speculative way. They used optical spectral data 
to probe the slowest ejecta in supernova explosions and em- 
ployed the radio observation to trace the fastest component 
of the outflow. They then took these two data points to es- 
timate the energy- velocity distribution. Their results may 
be biased because the fast moving material identified by ra- 
dio observation might be the decelerated GRB/XRF ejecta 
rather than the fastest component of the main SN explosion. 
If so, it is not a continuous distribution of matter between 
the two data points (Soderberg et al. 2006; Xu et al. 2008). 
Fairly speaking, observationally so far we do not have 
a reliable estimate of the initial kinetic energy- velocity dis- 
tribution of (most) HN outflows in the velocity range of 
fi ~ 0.1 — 0.5. Theoretically the standard hydrodynamic 
collapse of a massive star (Tan, Matzner & Mckee 2001) re- 
sults in a kinetic energy profile of the SN explosion E{> 



r/3) (X (r/3)- 



Such a steep function, however, is incon- 



sistent with the constraint from the optical data of SN 
20031w (Mazzali et al. 2006), for which a rough estimate 
gives E{> r/3) oc (r/3)"^'^. Motivated by this fact, we as- 
sume that all HNe associated with GRBs have a fairly steep 
energy distribution, which is generally written as i5(> r/3) = 
A[r/3/(r„/3o)]"'' for (3 < 0.5, where Fo = (1 - /3^)-i/^ For 
SN1998bw/GRB980425, SN2003dh/GRB030329, 



SN20031w/GRB031203, and SN2006aj/GRB060218, optical 
modelling suggests A ^ 0.2 — 6 x lO^^erg and VoPo ~ 
0.04 - 0.1 (Soderberg et al. 2006). The parameter q, how- 
ever, is not reliably determined in most cases. For sim- 
plicity, we approximate E{> r/3) = A\rf3/{ToPo)]^'' as 
E{> 13) = AiPlPo)-" for F ~ 1. 

In a SN explosion, the larger the velocity, the outer 
the ejecta. When the fast component is decelerated by the 
medium, the slower part will catch up with the decelerat- 
ing shock front. As a result, the total kinetic energy of the 
shocked medium increases and the deceleration of the shock 
is suppressed. In the quasi-similar evolution phase of the HN 
shock, the fastest component has swept enough medium and 
has got decelerated. Significant part of the initial kinetic en- 
ergy of the HN material E{> (3) has been used to accelerate 
the medium to a velocity ~ /3. So when we talk about the CR 
acceleration in the blast wave, the E{> (3) mentioned there 
actually represents the total kinetic energy of the shocked 
medium moving with a velocity /3. For a medium taking the 
profile n ex i?^*^ (0 < fe < 3), the rest mass swept by the HN 
blast wave is Mmed = /„ ^imrripR^dR cc R^~''. For l3 > /3c, 
conservation of energy gives E{> (3) ~ Mmod/3^/2, i.e.. 



P 



-(9 + 2) 



(xR' 



(1) 



With the relation that R ~ (31, the dynamics of the HN 
outflow is described by (for [3 > l3o) 



3-fc 
(3 (X t 5 + q-k ^ 



(2) 



In the next section we'll show how q and k influence the CR 
spectrum. 



2.2 The medium the hypernova outflow expands 
into 

As shown in Eq. (2), the dynamics of a HN shock de- 
pends on the medium profile sensitively. Here we review 
the medium profiles of all four GRB-HN events, based 
on the GRB and/or HN afterglow modeling. For GRB 
980425, the medium is Wind like (i.e., k ~ 2). The af- 
terglow modeling favors an unusual small A^, ^ 0.01 — 
0.04 (Li & Chevalier 1999; Waxman 2004), where A, = 
(M/1O"^M0 yr"^)(vw/10® cm s"^)"\ M is the mass loss 
rate of the progenitor, and v„ is the velocity of the stellar 
wind. For GRB 030329, the circumburst medium is found to 
be homogeneous (i.e.. A: ~ 0), as shown in many independent 
investigations (Frail et al. 2005; Pihlstrom et al. 2007; van 
der Horst et al. 2008; Xue et al. 2008). For GRB 031203, 
after modeling the radio data, Soderberg et al. (2004) got a 
constant ^ n ~ 0.6 cm"^. For GRB 060218, the high qual- 
ity radio data supports the homogeneous medium model 



^ There is an additional evidence for such a conclusion. With an 
electron energy distribution index p ~ 2.6 inferred from the radio 
data, the X-ray spectrum Fi, oc j/""-* (Soderberg et al. 2004) 
suggests that the X-ray emission is below the cooling frequency 
of the forward shock. For a free-wind medium, the late-time (t > 
1 day) X-ray light curve should drop with the time as t^^'^, 
deviating from the detected t~^ decline (e.g., Ramirez-ruiz et al. 
2005) significantly. If the medium is ISM-like, the expected late- 
time decline is t^^'^, consistent with the data. In this scenario, 
the early time (i < 1 day) X-ray flattening, as those well detected 



with n ~ 100 cm"^ (Fan, Piran & Xu 2006; Soderberg et 
al. 2006). As such, we have no compelhng evidence for a 
Wind-like medium surrounding most GRBs, even for those 
associated with HNe. The physical reason is not clear, yet. 
A post-common envelope binary merger model (e.g.. Fryer, 
Rockefeller & Young 2006) or a fast motion of the Wolf- 
Rayet star relative to the ISM (van Marie et al. 2006) may 
be able to solve this puzzle. 

Actually, a free wind medium, supposed to surround the 
progenitor, is unlikely to be able to keep such a profile up 
to the radius 



RicciPo) ~ 3 X 10'' ( 



lOMf 



-)A;^ 



(3) 



where Mcj is the rest mass of the GRB-associated HN ejecta, 
and 7?dGc is the deceleration radius. This is because during 
their evolution, massive stars lose a major fraction of their 
mass in the form of a stellar wind. The interaction between 
this stellar wind and the surrounding interstellar medium 
creates a circumstellar bubble (e.g., Wijers 2001; Ramirez- 
Ruiz et al. 2001; Dai & Wu 2003; Chevaher et al. 2004; van 
Marie et al. 2006). The analytical calculation suggests that 
the free wind of a Wolf-Rayet star usually terminates at 
(Chevaher 2007) 



Rt = 5.7 X 10^** (- 



-i-'mos 



p/k 



^-1/2 .1/2 



KU cm,(4) 



103 kms-i"l05cm3 K' 

where p is the pressure in the shocked wind and k is the 
Boltzmann constant. This is confirmed by observations of 
Wolf-Rayet nebulae, such as NGC 6888 and ROW 58, which 
also have radii of the order of a few pc (Gruendl et al. 2000) . 
Here we take the numerical example given in Figure 1 of 
Chevalier et al. (2004) to show that the HN outfiow is mainly 
decelerated in the ISM-like medium region. In their numeri- 
cal example, n ~ 0.5-R^g cm~^ for i?i8 < 1.2. The total mass 
of the free wind medium is thus ~ 7 x 10^"^ Mq ^ Mcj . Here 
and throughout this work, the convention Qx = Q/IO^ has 
been adopted in cgs units. 

We therefore conclude that the medium is most likely 
to be ISM-like at the radius where the HN outflow has been 
decelerated to /3 < 4/3o. This can be also understood as 
follows. One can infer from eq.(l) that for k — 2 the out- 
fiow component with a, (3 > l3o will decelerate at a radius 
RdeciP) « (/3//3o)-^'+'''7?dec(/?o). So for (5 ~ 4/3o and g ~ 2, 

-Rdec(4/3<,) ~ 4 X 10"^7?dec(/3o) ~ Rt. 



3 SPECTRUM OF COSMIC RAY PROTONS: 
OBSERVATION AND INTERPRETATION 

3.1 A new CR proton component in the energy 
range of 10^^ - lO^*"''' eV 

The spectrum of protons steepens suddenly at the first knee 
by a factor of 

A7(I) 2.1. 

In view that the spectra of heavier particles would steepen 
at higher energies, the likely interpretation of the steepening 



in Swift GRBs, can be attributed to an energy injection from the 
central engine. 



of all CRs at the first knee is the sudden decline of the light 
particles such as H and He (see Hillas 2005; Horandel 2008 
and the references therein). 

The proton CR spectrum before and after the second 
knee, after subtracting the modeled "galactic" component, 
can be roughly estimated as (Ulrich et al. 2004; Antoni et 
al. 2005; Hillas 2005; Horandel 2008) 



dEcB 



ex 



p-2.4 



for 0.1 < EcR.n < 3, 
for 3 < £cR,i7 < 30, 



(5) 



which indicates the factor of spectral steepening is A7(II) ~ 
—0.9. Above the ankle, the CR spectrum changes to S^p^' , 
so the factor of fiattening is A7(III) ~ 0.6. 

The interpretations of spectral changes at the second 
knee and at the ankle are much less clear. Hillas (2005) in- 
terpreted them as a result of an extragalactic component 
with a spectrum oc E^^"^ suffering losses by the interaction 
between cosmological microwave background radiation and 
starlight. In this work, we consider the detected spectral 
change around the second knee is due to the energy- velocity 
distribution of HN outfiows. 



3.2 Theoretical interpretation 

3.2.1 Self-magnetized shock acceleration model 

In this model, the magnetic field of the upstream is assumed 
to be amplified significantly by the CRs themselves (e.g., 
BeU & Lucek 2001). 

With Eq. (2), we have the radius of the forward shock 

2 + 9 

front as _R ex t5+<j-fc The maximum energy accelerated by 
the forward shock can be estimated by (Bell & Lucek 2001; 
Berezhko & Volk 2004; Ptuskin & Zirakashvili 2005) 

fc-4 + g(2-fc)/2 4-fe-q(2-fe)/2 

Smax(fc,g)~^/3eBi?oci 5+F^^ oc/3 ^^^ , (6) 

where B ex /Si?^*'' is the magnetic field in the upstream of 
the shock, which is on the same order of that of the shocked 
medium. 

In the ISM case (i.e., fc = 0), we have 



i5max(0,g)oC/3(^-«>/3^ 

while in the WIND case (i.e., k — 2), 

S,nax(2,q)cx/3'. 



(7) 



(8) 



Here we do not present the numerical coefficient of 
Eyna.^{k,q) because Wang et al. (2007) and Budnik et al. 
(2008) have already shown that for typical parameters, a 
/3 ~ /3o ~ 0.1 is high enough to accelerate protons up to 
~ 10^^ eV regardless of k. In a stellar wind medium, the 
HN shock front with a /3 ~ 4/3o can accelerate protons up to 
~ 3 X 10^* eV. Therefore the CR protons above the second 
knee but below the ankle are mainly accelerated by the HN 
shock in the velocity range of ~ /3o — 4/3o. 

To get an estimate of the spectrum of the accelerated 
particle, following Berezhko & Volk (2004) and Ptuskin & 
Zirakashvili (2005) we assume: (1) The particles with an 
energy iSmax escape the shock immediately; (2) The total 
energy of the accelerated particles at an energy i?cR = 
i?max(/3) is proportional to E{> 0). In view of the rela- 
tions E{> 13) (X [£max(/3)]"*'"'<''"''' for fc = and E{> f3) ex 
[Smax(/3)]"'''" for fc = 2, we have 



dN 



dE, 



CR 



CR 

-(2+<S)-f 
CR 



for A; = 0, 
for A: = 2, 



(9) 



where S ~ 0.4 is introduced to account for the proton spec- 
trum in the energy range of 10'^'' — 3 x 10^^ eV. As /3 < /3o, 
E{> (3) (X 13 if the energy loss of the HN shock is ignorable. 
The accelerated proton spectrum should be oc E^^ . This 
answers why there comes a spectrum change around the sec- 
ond knee if -Emax(/3o) ~ 3 x 10^'^ eV. 

With a 5 = 0, to match the detected proton spectrum 
dN/dEcB. oc Ef^^'^, one has to have q ~ 2.6, which is very 
close to that of SN 20031w and SN 1998bw reported in Soder- 
berg et al^ (2006). Therefore Wang et al. (2007) concluded 
that the self-magnetized shock acceleration model could ac- 
count for the spectrum data. However, a few puzzles have to 
be solved before accepting this argument: (I) If 5 = 0, some 
novel effects are needed to interpret why the proton spec- 
trum departs from E^^ significantly in the 10^® — 3 x 10^^ eV 
range. The authors also need to explain why these effects, if 



any, disappeared in the 3 x (10 



10^ 



eV range. (II) A 



wind profile holding to a radius ~ 1 — 10 kpc is crucial for 
their argument. If the medium is ISM-like when the outflow 
gets decelerated to /3 < 0.4, Wang et al. (2007) 's approach 
would yield a spectrum 

sfr^'^cJ (10) 

for g ~ 2, which is too steep to be consistent with the data. 
We take this puzzle as a potential limit of their interpreta- 
tion. 

Let's investigate whether a specific wind-bubble can 
solve this puzzle. We assume that the free wind profile is 
terminated at a radius ^ Rt ^ 10^^ cm and is followed 
by an ISM-like shell. Suppose that the shell is so mas- 
sive that the deceleration of the whole HN outflow occurs 
at i? ~ i?t ~ const, we have iSmax ~ Z[3eBR oc [3B. 
If the shell is not dense enough to form a strong reverse 
shock, i.e., the forward shock velocity decreases continually 
rather than abruptly, then B oc /3n^". As a result, we have 



P'^n^/'^ and dN/dE^^ oc E, 



-(2+<s)-f 

CR ' 



provided that 



the CR protons in the energy range of ~ 3 x (10 , 10 ) 
eV are mainly accelerated in the shocked shell. Though 
such a possibility is attractive, the request that the re- 
verse shock does not form is hard to satisfy. This is be- 
cause at a radius ~ i?t ~ 10^® cm the number density of 
the wind medium riw ~ 3 x 10^* cm^'* A,,_i_R^j^g. On the 
other hand, the assumption that 47ri?t ntnip ~ Mcj requires 
that nt ~ 1 cm~^ (Mcj/lOM0)i?~ig. So we have a den- 
sity contrast nt/ww ~ lO"'. The forward shock expanding 
into the dense shell will have a pressure ~ P^ntm^c? /Z, 
which is much higher than that of the shocked wind medium 
(~ (i^nv,m-p<? /?>). A pressure balance will be established 
by a strong reverse shock penetrating into the shocked 
wind medium. Therefore the forward shock velocity is much 
smaller than /3 | shocked wind medium aud can not accelerate 
protons to an energy ~ 10^* eV. The reverse shock with a 
velocity f3r ~ l3 | shocked wind medium plausibly plays a more 
important role in accelerating high energy CR protons. The 



* Please see section 2.1 for the discussion of uncertainty of the q 
obtained in their way. 



shocked wind medium has only a very small mass (relative 
to Mcj). The reverse shock gets weak after penetrating into 
the dense HN outflow which has a density comparable to rit. 
Then the forward shock velocity increases and significant CR 
acceleration in the forward shock front is possible. A detailed 
numerical calculation, like Ptuskin & Zirakashvili's (2005), 
is needed to draw further conclusions. 



3.2.2 Gyroresonant stochastic acceleration model 

The maximum energy-gain rate due to the stochastic Fermi 
acceleration for marginally relativistic shock can be esti- 
mated as (Dermer 2001b) 



dE^ ^ etn.b(^-l) . 

dR 2»/2 ^ ^ 



,1/2 



Ecn 



ZcB^JaRP 



(11) 



where Z is the atomic number, Sturb is the ratio of plasma 
turbulence to the shock energy density, _B* « 0.4 n^'^ej 
Gauss, /a ~ 1/12 is the ratio of the width of the swept 
medium by the shock to R (Dermer & Humi 2001), and v 
is the spectrum index of the turbulence {v — 5/3 for Kol- 
mogorov turbulence and 3/2 for Kraichnan turbulence). 

Dermer (2001b) took a /3 ~ const., integrated eq.(ll) 
over R, then got £'max(fi). However, currently P evolves with 
R. As shown below, the smaller the radius, the larger the P 
and the higher the Bmax. Very energetic CRs can be accel- 
erated at early times but can not be accelerated continually 
because of the adiabatic cooling. Taking into account the 
adiabatic cooling effect, eq.(ll) takes the new form 



O'-E^CR ^ eturb(t'- 1) ^^f, ^2^ 2^'^£;CR -^v-1 Sen ^2) 



dR 



23/2 



-ZeB^P^i 



ZeB,fARP 



R 



Now i5max can be estimated by setting ^^"- = 0, then we 
have 



P ^ r £turb(i' - l)P ^i/(2~v) ZeB^f^RP 

2/A V2 

ISM medium. In this case, we have 

1/2 1/2 



(13) 



Smax(ISM) ~ Zn'o'€'j^;'_^R 



10« eV (^ypt,{l2U)-^ {v = 5/3), 
10" eV {^fp^_,{12U)-^ {v = 3/2). 



(14) 



The energy conservation 4,-k R^ P^ nrUpC? / 'i ~ E{> P) 
yields R « IQ^"" cm A^/'^/SZ^'^+'^/'n, 
Eq. (14), we have 



19 „„ 4I/3 «-(.+2)/3„-l/3 Cojnbij^ijjg „ith 



£max(ISM)0C/3' 3(2-*') "\ (I5) 

i.e., £;max(ISM) oc /jC^-?)/^ for v = 3/2 and oc /3(i«-9)/» for 
V = 5/3, both are sensitive to /3. 

WIND medium. In the termination wind shock model, 
the stellar wind profile may hold up to a distance ~ 10^* 
cm (e.g.. Chevalier et al. 2004). In this case, n = 3 x 
10^^ A-^R-^ cm^^ Now B, ^ Q.2A\{\e](\R-^ Gauss and 



-Emax(wind) ~ ZA 



1/2 1/2 



2 X 10^5 eV (^)3/31i(12/a) 

2 X 10« eV [^-f^fp^_,[l2fA)-' {v = 3/2) 



-; irJillm 



As shown in Eqs. (14) and (16), for eturb ~ 0.5 and v — 
(3/2, 5/3), at /3 ~ /3o ~ 0.1, we have S^^x ~ (10^^, lO^'')^ 
eV. Below we focus on the case oi v — 3/2, because in the 
caseofv = 5/3 the request of i5niax(/3o) ~ 3x10^^ eV is more 
difficult to satisfy. For /3 ~ 0.5, the stochastic gyroresonant 
acceleration is able to accelerate protons to ~ 10^^ eV (see 
also Dermer 2001a). The accelerated particle spectrum is 
thus (« = 3/2) 

dEr°' 



E^CR 



-^CR 
-^CR 



7-g 



for fc = 0, 
for fc = 2. 



(17) 



As shown in section 2.1, the main deceleration of the HN 
outflow is very likely to be in an homogenous medium. 
The accelerated protons have a spectrum dN/dE^^ ex 



E, 



-(2.4+3q/(7-q)) 
CR 
3q/{7 - q) 



To match the observation A7(II) 
-0.9, we need 



g~ 1.6. 

This is surprisingly close to the value --^1.7 that is inferred 
from the optical modeling of SN 20031w. Detailed optical 
modeling of more GRB-associated HN explosions is highly 
needed to better constrain q and then confirm or rule out 
our interpretation. 

If GRB-associated HNe expand into a Wind bubble- 
like medium, a flatter CR spectrum in the higher energy 
range would appear. At a small radius (say, < 10^* cm), 
the medium is free Wind-like and the accelerated particle 
spectrum is ex i?cR ' which then gets steepened by a 

factor of q{2 + q)/[3{7 — q)] ~ 0.4 for q ~ 1.6 after entering 
the ISM-like medium. Such a flattening seems not enough 
to match the observation A7(III) ~ 0.6. So CRs above the 
ankle may be mainly from AGNs, as indicated by the recent 
analysis of the correlation of the highest-energy CRs with 
nearby extragalactic objects by the Pierre Auger Collabora- 
tion (Abraham et al. 2007). 

The rate of local GRB-associated HNe only accounts for 
~ (0.1 - 0.5)% of that for aU local SNe (Delia Valle 2006; 
Soderberg 2007). The typical energy of these HNe, however, 
is tens times larger than that of the normal SNe. Roughly, we 
expect that a fraction ^ 10% of CR protons at 3 PeV could 
be attributed to GRB-associated HNe. It is enough to match 
the observation (Ulrich et al. 2004; Horandel 2008). So the 
CR proton spectrum in the energy range of 10^® — 10^*'^ eV 
may be quantitatively interpreted. 



4 DISCUSSION AND SUMMARY 

The particle acceleration in marginally relativistic HN 
shocks are discussed. The GRB-associated HN outflows are 
assumed to have a fairly steep energy distribution against 
their velocities, i.e., E{> /3) oc /S"' for q ~ 1.7, as in- 
ferred from the optical modelling of SN 20031w (see sec- 
tion 2.1 for details). A significant fraction of a HN's kinetic 
energy is carried by the material moving with a velocity 
> l3o{~ 0.1), driving an energetic shock wave into the sur- 
rounding medium. The cosmic ray protons above the sec- 
ond knee but below the ankle may be accelerated by the 
HN shocks in the velocity range oi P ^ (1 — 4)/3o. To satisfy 
this velocity bound, the HN outflows associated with GRBs 



must have reached a very large radius where the surround- 
ing medium is very likely to be ISM-like (see section 2.2 for 
details). With this argument, the self-magnetized shock ac- 
celeration model adopted in Wang et al. (2007) would yield 
a very steep spectrum that is inconsistent with the observa- 
tion unless the medium surrounding the HN is a free wind 
holding up to a radius Rdec ~ 10 kpc (Moj/IOMq)^^^^. 
Such a request seems difficult to satisfy. A highly specula- 
tive solution is that the particle acceleration mainly occurs 
in a narrow "dense" shell that terminates the free wind at 
a radius ~ 10^^ cm (see the last paragraph of section 3.2.1 
for details). 

In this work, we find that for q ~ 1.6, the stochastic 
gyroresonant acceleration model can account for the spec- 
trum change of high energy protons around the second knee 
(see section 3.2.2 for details). As a consequence, the stochas- 
tic gyroresonant acceleration mechanism in relativistic GRB 
forward shock may account for part of the ultra-high energy 
CRs (~ 10^° eV), as suggested in Dermer (2001b, 2007) and 
Dermer & Humi (2001). A typical 



-7A7(II) 
- A7(II) 



1.6 for A7(II) 0.9, 



if confirmed in future optical modelling of the GRB- 
associated HN explosions, will be a crucial evidence for our 
current speculation. 
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